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Abstract 

We consider the possibility that neutrinos are coupled very weakly 
to an extremely light scalar boson. We first analyze the simple problem 
of one generation of neutrino and show that, for ranges of parameters 
that are allowed by existing data, such a system can have serious conse- 
quences for the evolution of stars and could impact precision laboratory 
measurements. We discuss the extension to more generations and show 
that the general conclusion remains viable. Finally, we note that, should 
such a scalar field be present, experiments give information about effective 
masses, not the masses that arise in unified field theories. 



1 Introduction 



In the last few years, it has been suggested that neutrinos might intcract weakly 
among themselves through the exchange of a very light scalar particle jl], ^ , with 
possible consequences for the evolution of the Universe and for the propagation 
of neutrinos from distant events. In many of these discussions, one assumes that 
ncutrinos are distributed according to the usual Big Bang scenario and Standard 
Model physics, the effects of scalar exchange being treated as a perturbation. 
In this paper we examine that assumption. 

This problem is a special case of the general problem of relativistic fermions 
interacting through the exchange of scalar and vector bosons, and the gen- 
eral formalism has been worked out under the name Quantum Hadrodynamics 
(QHD) and applied to the study of nuclear physics ||. We shall show that, for 
a wide range of parameters, neutrinos will tend to cluster, that these clustcrs 
could be of a size to affect stellar formation and dynamics and that there may 
be observable consequences for physics within the solar system. In fact, if the 
clustering is strong enough so that the density of ncutrinos within a cloud is 
sufficicntly large, terrestrial laboratory experiments can be affected. 

Over the last decade, many groups studying the endpoint of the Tritium 
beta ray spectrum for signs of neutrino mass have reported a best fit value for 
the square of the anti-neutrino mass less than zero § |, |, @, |, §]. Robertson 
et.al Q point out that this result could be obtained by assuming a very small 
branch due to the absorption of relic neutrinos, provided the density of such 
neutrinos were some 10 13 higher than usual cosmological values. In choosing 
parameter ranges for the examples in this paper, we have kept this idea in mind. 
We must emphasize, however, that the general feature of cloud formation and 
its consequence for the evolution of structure in the history of the Universe is 
quite robust, and must be considered whatever the eventual resolution of the 
anomaly in Tritium beta decay. 

In this work, we only consider the effects of light scalar exchange. The effects 
of the known vector exchange (Zq) are far too small (really, too short ranged) 
to affect these results. The exchange of a light vector particle is severely con- 
strained by data. To avoid problems with the axial anomaly (the neutrinos do 
couple to the Zq) one must either invent many new fermions or demand that 
the light boson couple to known leptons or quarks, which quickly leads to con- 
flict with experiments designed to test for a fifth force jn|. Furthermore, the 
self shielding of a vector exchange, while it might allow for the development 
of a ncutrino-antincutrino plasma, will not allow for the coherent action re- 
quired to drive cloud formation. Thus, we only treat scalar exchange. Even so, 
there remains the rich possibilities of different couplings to different generations. 
Aside from a few comments, we leave that to further work, concentrating our 
discussion on the simpler system of one flavor of neutrino. 

The paper is organized as follows. In section 2 we review the treatment of 
infinite matter in QHD and apply that to the problem at hand. In section 3 
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we do the same for finite clouds of neutrinos. In section 4, we describe the 
consequences such clouds would have on the the evolution of structures in thc 
Universe and, in section 5, we confront this picture with what data exists, 
extracting limits on the parameters. In section 6 we discuss the changes in the 
analysis of Tritium beta decay experiments within such neutrino clouds, and 
comment on the effects of such clouds on stellar dynamics in section 7. In section 
8 we discuss some aspects of the extension to include more than one generation, 
illustrating the remarks with special cases applied to two generations. We offer 
our conclusions in section 9. 



2 The Infinite Problem 

The Lagrangian for a Dirac Field interacting with a scalar field is well known: 

£ = $(ip - m v )i + - [cj)(d 2 - m 2 s )(f\ + g^cj) (1) 

which gives as the equations of motion 

[5 2 + m 2 ]^ = g^ (2) 

[ip - m v \ i' = -g<H>. (3) 

We omit nonlinear scalar selfcouplings here, even though they are required to 
exist by field theoretic selfconsistency, as they may consistently be assumed to 
be sufficiently weak as to be totally irrelevant. 

We look for solutions of these equations in infinite matter which are static 
and translationally invariant. Equation (^) then gives 

= -^, (4) 
mj 

which, when substituted in (^|) gives an effective mass for the neutrino of 



m,. = m 



2 



:H>- (5) 



where m v is the renormalized vacuum mass that the neutrino would have in 
the absence of other physical neutrinos. 

These equations are simply the equations of Quantum Hadrodynamics J3|, 
and we will be using them in a small coupling regime where there is no question 
of the validity of neglecting higher order processes. 

These equations are operator equations. We next act with each of these 
equations on a state |fi) defined as a filled Fermi sea of neutrinos, with a number 
density p, and Fermi momentum kp, related as usual by p = kp 3 /(67r 2 ). Thc 
operator ipip acting on this state gives 



^|fi) = ^/ d 3 k =| |fi), (6) 
( 27r ) 3 J\k\<k F y/(m* v ) 2 + k 2 
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where w is the number of neutrino states which contributc — w = 2 for Ma- 
jorana neutrinos and w = 4 for Dirac neutrinos. Thus the effective mass is 
determined from the integral equation 

(o) 9 2 w f kp 2 m* 



To discuss the solutions of this equation we reduce it to dimcnsionless form, 
dividing by m v °\ and introducing the parameter K = 9 ^T^i ^ , and thc 
variables y = -^j,x = —jgj,x F = -^foj- Then equation (Q) becomes 



i-f^o/ 7= (8) 
•'o \/y + x* 



yK 



e fx f — y hi 



e F + x F 
V 



(9) 



with e F = \J x f + V 2 - One can regard equation (j^) as a non-linear equation 
for y as a function of either e F ot x F . As a function of ep, y is multiple valued 
(when a solution exists at all), whereas y is a single valued function of :tf- 

The total energy of the system is a sum of the energy of the neutrinos, 
E v = e„ mi^ wN, and the energy in the scalar held, E s = e s m v wN, where 
N is the total number of neutrinos in each contributing state. These expressions 
serve to define the per neutrino quantities e v and e s . Also, E s = £ s V , where 
£ s = \m 2 s 4> 2 is the energy density of the (here uniform) scalar field. 

One finds that 

e v = — / x 2 dx\J x 2 + y 2 
X F Jo 

3 f x F e F ^ x F y 2 e F y 4 ^ / e F + x F \ \ ^ 



and 



K 3 2 / f XF x 2 dx 
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3-y 

1 3 



2 X F \Jo \J x2 + y 
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2K x 3 F 

Notice that for large values of x F 

y -> 



(i-j/r- (ii) 



K x% 

3e F 

~4~ 



3 



3xf 
3 1 

Cs ^ ~~ 
It is also useful to note that, for small Xp, 

K x 



(12) 
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es - f f . (13) 

For the neutrino system to be bound, the minimum of e = e„ + e s as 
a function of density (or xf) must be less than 1, its value in the zero density 
limit. 

These points are illustrated in Figures 1 and 2, for the case Kq = 20. In 
Figure 1, the average energy per neutrino, e, and the Fermi energy, ep, &re 
plotted, in units of the vacuum mass mo, against thc Fcrmi momentum in the 
same units, xp. In Figure 2, the same information is displayed as e vs ep, 
demonstrating the multivaluedness of the solution viewed this way. The general 
shape of Figure 1 will be exhibited for any value of Kq, but only if Kq is large 
enough will there be a minimum with e < 1. This is displayed in Figure 3, from 
which we may deduce that Kq must be greater than w 3.3 for binding to occur. 
Note also that the local minimum disappears for Ko = 2.67, showing that 
there is no metastable state for smaller values of the effective coupling strength. 



3 Neutrino Clouds 

Note from equation (|l3|) that the energy varies as 1 + const. x kp near xp = 0, 
where const. is positive. Even with K large enough to produce binding, at low 
density e will be greater than 1, then will decrease with increasing xp to its 
minimum value before increasing as xf — * oo. Thus, for low densities, it may 
be energetically favourable for the neutrinos to form finite size clouds. 

To investigate this possibility we use the Thomas-Fermi approximation to 
obtain the scalar ficld through thc equation 

V 2 ^ + m 2 J = g(^), (14) 

substituting the local density value for ('ipip). 
Noting that 

y = l -^m^ ( 15 ) 
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we convert the local density equation ( |l4| ) for (j> to a local density equation for 

G{y,e F ,K ) (16) 



d 2 y 2 dy 

dz 2 z dz 



with 



X F 



G{y,e F ,K ) = -l + y[l + K , (17) 

and z = r m s , where r is the radial distance from the center of the cloud 
and we have assumed spherical symmetry. In this equation, e F (which is a 
chemical potential for thc ncutrinos) is constant with position, but y and thus 
x F = \J e p ~ y 1 are functions of position. Clearly the surface of the cloud is 
defincd by the condition that xf = 0, as there the neutrino density drops to 
zero, at the (scaled) radius z = z . For z > z , x F remains zero, although thc 
scalar field does not immediately vanish. Instead, outside the cloud the solution 
of the differential equation for the effective mass y becomes 

g — (»— zo) 

y =l-{l-e F )—j— (18) 

z/z 

which satisfies the condition y{z ) — e F at the surface, and the condition that 
y — > 1 as z — > oo. From this exterior solution we compute the values of y and 
at z = z , to which the interior solution must be matched: 



y{z ) = e F 

djl 
dz 



In the usual way one solves the interior equation numerically starting from an 
assumed central value, y(0), which is adjusted to satisfy the matching conditions. 
While this technique works in a straightforward way for small clouds, as the 
cloud increases in size one needs to know y{0) to greater and greater precision, 
and needs to adopt a special technique to handle such cases numerically. 

For large values of z the central region will approximate the conditions of 
the infinite system, and so we expand the function G of equation (|l^) as a Taylor 
series in y about y , the solution of G{y,e F , K ) = 0, which is the effectivc 
mass in the infinite system. Explicitly 

G{y,e F ,K ) 

l yo 

K e F ( y x \ [ 
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3- \{y-y y + .... (20) 
a^o yoj } 

Keeping just the first term in equation (pfj|), thc equation (|l^) has the solu- 
tion 

. sinh(Kz) 

y~y = A 2 1 

nz 
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where 

n 2 = J-2 + — -K e F x \. (22) 

(It is straightforward to show that k 2 > for the cases of interest, which are 
those for which e is close to the lower branch in Figure 2.) 

The numerical solution of equation ( fL6| ) is started at z = z\ > 0, with an 
assumed value y(zx) = y + Ay, where Ay is small enough that the higher 
order terms in equation (^0|) can be safely neglected. The value of Z\ is adjusted 
to satisfy the outer matching conditions. One can then use equation (^lj) to 
extrapolate y(z) from z = z\ to z = 0, obtaining 

2/(0) = yo + . Ay. (23) 

smn(KZi) 

In this way we can set y(0) at the extraordinary level of precision required to 
satisfy the matching conditions. 

The general form of the solutions does not depend on the particular value of 
K . An example of solutions with different numbers of neutrinos is presented 
in Figure 4 for K = 20. To construct the Figure, values of e F were chosen 
and a solution obtained, giving the density as a function of the radial distance 
from the center of the cloud. That density was then integrated to obtain the 
total number N displayed with the appropriate curves. The units are such that 
N=l would correspond to (mJJ/m s ) 3 total neutrinos. For the Thomas-Fermi 
approximation to work well, that number should be large compared to 1, which, 
as will be seen below, is the case. For later reference, we display the radial form 
of the gravitational force, f g , due to such a distribution in Figure 5. This force 
is presented in arbitrary units, since the masses have been scaled out, and is 
calculated as f g = 72 J e z 2 dz. Since, for finite clouds, the scalar field is 
no longer always constant, there is an additional contribution|| to the scalar 
field energy density, i(V</>) 2 , which has been included in our calculation. 

To understand the distribution of cloud sizes to be expected, one needs to 
know the average energy per neutrino as a function of size or of total number. 
For K = 200, Figure 6 displays both the logarithm of the cloud radius (i.e. 
the radius at which the density goes to zero) and the logarithm of the difference 
between the average energy in the cloud and the average energy in infinitc 
neutrino matter at its minimum for the same value of K versus the logarithm 
of the total number. The points are the results of calculations corresponding 
to the indicated values of e F . The line through the radius values is a straight 
line with slope |, that through the energy differences has slope Evidently, 
for large N, the radius scales as N 1 ^ 3 and the energy difference scales as N^ 1 / 3 , 
appropriate to a surface tension, indicating a preference for large clouds. Thc 
line through the e F values is to guide the eye. 
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4 Early Universe 



Consider the effccts of such clustering on the evolution of structures in the 
early Universe. Throughout the following discussion, we assume that K is 
large enough to produce bound systems. At an early enough epoch the density 
will be sufficiently high that the effective mass is negligible. At that epoch, 
there is no difference between the interacting neutrinos and the relativistic, 
non-interacting gas that is usually assumed. Conscqucntly, these ncutrinos will 
cxpand and decrease in density according to the standard scenario until the 
increase in the effective mass begins to make a difference, which will occur at 
about the value of the density (xf) corresponding to the minimum energy per 
particle for infinite matter. 

Were it possible to remove energy from thc ncutrino gas and entirely from the 
Universe, so that the gas could be viewed as having zero temperature, that would 
be the end of the discussion. We would, however, be left with a conundrum. Thc 
ncutrinos could tolcrate no furthcr expansion but the Universe, being driven by 
all sources of energy density, would be required to continue to expand. This 
would result in one neutrino cloud located somewhere (defining a "center" ?). 

That, of course, is not the situation. The ncutrinos will have a temperature 
comparable to that obtained for an expanding, non-interacting gas. As the 
expansion continues, that temperature will be converted into (effective) mass 
and the gas will become supercooled, followed by fragmentation into clouds. 
Note that no additional dissipation is required, unlike the case where clouds 
coagulate from free particles. The point here is that the neutrinos were born 
within a cloud and never achieve a state in which the effective mass rises to its 
vacuum value. 

Many ncutrinos have been born at later times through normal stellar burn- 
ing, supernova explosions or other processes. When they encounter a neutrino 
cloud, the coherent forward scattering is easily large enough, even for very small 
values of the coupling to the scalar field so that individual scatterings are small, 
to cause the neutrinos to lose energy through the Bremsstrahlung of scalars, 
providing additional dissipation. 

Two factors drive the size distribution of these clouds. The first is thc 
distribution of fluctuations, which we assume follows Harrison-Zel'dovich |p"l|| . 
The second is the increase in energy per neutrino with decreasing cloud size 
discussed in section 3. The latter effect provides for a mechanism to cut off 
the distribution of cloud sizes below some smallest value, the actual efficacy 
of which depends on the detailed parameter values. The general form of thc 
distribution is 

P(N) oc N~ 2 exp(-C/iV 1/3 ). 

Should this process occur before recombination, which, as we shall see in 
Section 5, is not beyond reason, then the existence of neutrino clouds would 
have a profound effect on the evolution of small size structures. (By small, 
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in this context, we refer to structures of the size of solar systems, stars or a 
bit smaller.) At recombination, when matter decouples from the background 
photon gas, there will be a pre-existing collection of gravitational sources. Thc 
longer the time between cloud formation and recombination, the more these 
will appear to be point sources, but that does not strongly affect the following 
argumcnt. Whatever the spectrum of fiuctuations in the baryon distribution, 
these pre-existing sources will nucleate baryon condensation with a distribution 
that more or less follows the size distribution of thc clouds. 

Many of these collections of baryons will be large enough to initiate nuclear 
burning and become stars; others will not. Of the latter, some will attract more 
baryonic matter from the ejecta of exploding stars to form later generation 
stars, while others will remain too small to evolve into stars and can provide 
cold, massive objects. Note that, even if a given cloud does not attract a 
compliment of baryonic matter, it will still function as a gravitational source. 
In either case, the increase in the energy per neutrino with decreasing cloud 
radius, discussed above, will provide a lower limit to the distribution of system 
sizes. Thus, the existence of neutrino clouds can serve as a seed mechanism for 
stars and could provide a similar seed for (or be themselves) smaller objects 
such as MACHOs@. 

This scenario suggests that all stars will have their associated neutrino cloud, 
not because stars attract neutrinos but, rather, because stars form within the 
gravitational well provided by pre-existing neutrino clouds. One may then ask 
if, during subsequent evolution, the star and its cloud remain together or if the 
star, buffeted by forces which ignore the neutrinos, is stripped away leaving 
the cloud to catalyze another object. This is a quantitative question which we 
discuss in section 7. 

5 Numerical Considerations 

Following the spirit outlined in the introduction, we shall consider here the 
constraints that may be placed on the parameters of the theory under the as- 
sumption of one surviving neutrino species. We shall take up the question of 
other generations in section 8. 

First, to impact terrestrial Tritium beta decay experiments, the density of 
electron neutrinos needs to be about 10 15 /cm 3 . Robertson et. al. |Q] reported 
that their data could be fit by assuming an additional branch, 10~ 9 of the 
normal decay, which fed each final state in the same proportion as the normal 
decay. Applying standard formulas |l3| for inverse beta decay, this rate requires 
a density of 6 x 10 15 /cm 3 . A recent paper by the Moscow group |j| indicates a 
preference for a 6x 10 -11 branch, which would imply a density of 3.6 x 10 1 /cm 3 . 

For the first case, 

k F = (6tt 2 p) 1/3 = UeV/c 
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while the second would imply 

kp = 5.5eV/c 

Throughout this section we shall use 10 eV/c, since we are arguing here for the 
scale of the parameters. (An actual fit requires a more complctc trcatmcnt of the 
modification of the spectrum, which we discuss further in section 6). This value 
of kp corresponds to a density of 2 x 10 15 /cm 3 , which is some thirteen orders of 
magnitude above that obtained in standard cosmology, assuming that neutrinos 
are uniformly distributed in space. If we assume that all neutrinos created when 
the weak interaction was sufficiently strong to keep thermal contact between 
neutrinos and photons survive to the present epoch, this implies that only 10~ 13 
of space is occupied by neutrinos. 

As a consequence, neutrinos propagating over galactic or inter-galactic dis- 
tances will, in the main, be free of the influence of condensed neutrinos and will 
propagate with their vacuum mass. Therefore, mass limits obtained from the 
spread of arrival times for neutrinos from SN1987a will apply to the vacuum 
mass jl4|. While lower limits have been obtained by some authors, we shall use 
a limit of 50 eV here. Coupled with kp = WeV, this implies 

x F > 0.2 

or 

K < 4000 . 

From the requirement that the system be bound, we obtain the lowcr limit 

K > 3.3 . 

For one generation, particle and antiparticlc (chiral left and right, if Majorana), 
Ko is related to the basic parameters by 

K = 45/tt^ 2 (24) 

where 

li = m s /mW (25) 

and 

a = g 2 /4ir . (26) 

To obtain limits on a directly, we consider those processes involving either 
the scattering of neutrinos through scalar exchange or processes involving the 
Bremsstrahlung of scalars. For the purpose of these estimates, we use only lead- 
ing terms, and it does not matter whether the neutrino is a Dirac or Majorana 
particle. Consequently, many existing limits on the coupling of a Majoron to 
neutrinos can be taken over directly. 
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In a study of thc limits on such a coupling in a number of elemcntary particle 
processes p5[ , it was found that thc most stringent limit was that given by 
Barger, Leung and Pakvasa |l6| from a study of (K12) decay. This gave 

a < 1CT 6 . 

A more stringent limit can be obtained from recent studies of nuclear doublc 
beta decay Jjj]] looking for the spectral signature which would accompany the 
emission of one extra particle. 

a < 1CT 9 . 

Note that this limit only applies if neutrinos are Majorana particles. 

For y/s 3> m v the total neutrino-neutrino scattering cross-section is given 

by 

er w 5na 2 /s . 

The strongest constraint from scattering comes from the neutrinos produccd 
by SN1987a. If we assume that all neutrinos born in the Big Bang survive 
and are distributed roughly proportionally to the baryons, then there are 10 66 
associated with the solar system. At 2 x 10 15 /cm 3 , that implies a radius of the 
cloud of » 5 x 10 16 cm. On this scale, the Earth is essentially at the center of 
the cloud. Thus the mean free path for 20 MeV neutrinos coming in from the 
supernova must be greater than 5 x 10 1 cm or the cross section must be less 
than 10 _34 cm 2 . For a mean neutrino energy of 7.5 eV (3kp/A) and a beam 
neutrino of 20 MeV, s w 1.5 X 10" 10 GeV 2 , giving a w 10" 17 cm 2 5 2 This then 
requires that 

a < 3 x 10~ 9 . 

Other processes, such as the observation of solar neutrinos or the observation 
of neutrinos from high energy accelerators lead to less stringent upper bounds. 

A more interesting bound comes from the consideration of the survival of Big 
Bang neutrinos to the present epoch. The dominant disappearance mechanism 
is for a neutrino and an antineutrino to annihilate into a pair of scalars. For 
Majorana neutrinos, this simply means that a chiral left and a chiral right 
neutrino transmute into a pair of scalars. The rate is 

uo w — a kp m . 

8 m s 

The total disappearance probability is the integral of this rate to the present 
from whatever initial time is appropriate. The time before neutrinos condcnse 
into clouds is irrelevant to this question as, throughout that time, the scalar 
field and the neutrinos are both expanding relativistic gases which were once in 
thermal contact. All that a large rate for this process can do is maintain that 
thermal contact since, under these conditions, the back rate is the same. Conse- 
quently, we consider the integral of the rate from a time when the temperature 
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Table 1: Infinite Matter Values 



A 
20 
200 
2000 



< e > 
.704 
.410 
.233 



Xp 

675 
405 
232 



y 

.2 
.06 
0186 



xf/v 



.3 
.15 
.08 



corresponds to about 10 eV, or t — 10 13 s. The bound is only weakly depen- 
dent on the assumption of radiation or matter dominance for the expansion. 
Requiring a decrease of less than 0.1 in ln(7V/iVo) yields 



It is tempting to consider scenarios with a above this limit but still below 
the limits imposed by scattering. In that case, one might argue that neutrinos 
are Dirac particles and that there is a fractional particle excess as for any other 
fcrmion. It is, in fact, such a scenario that would give only an excess of counts 
at the end points of beta decay spectra (see the discussion in the next section) . 
However, the limits on neutrino generations from nucleosynthesis calculations 
are based on the assumption of only left handed neutrinos and right handed anti- 
neutrinos. The scalar interaction considered here will equilibrate left and right 
handed neutrinos while they are in equilibrium through the weak interaction, 
leading to a factor of two too many states per generation. As the remarks above 
regarding the limits of the time integration for disappearance point out, it does 
not seem possible to try to argue for the disappearance of some states at an early 
cnough epoch that it could affect the neutron-proton equilibrium and, hence, 
nuclear abundanccs. 

At this juncture, one may question one's ability to accommodate the scalar 
field in light of the same argument. This question was addressed several years 
ago by Kolb, Turner and Walker [[l8). The results of that paper indicate that 
a very light scalar (the case here) would behave as an extra half a generation, 
which can be accommodated. 

Returning to the issue of Majorana neutrinos, the search for the evidence 
for their existence in nuclear double beta decay is a thriving industry |l9|] . 
If such decays were mediated by a mass term, modern experiments limit that 
mass to less than 1-2 eV. When taken with our standard value of Uf = 10 eV, 
this implies that u/xf < 0.2 independent of the vacuum mass. Table 5.1 lists, 
for some representative values of Kq, values of Xf, y and the ratio y/xF for 
infinite matter at its minimum energy per neutrino. This suggests a preference 
for larger values of K$ within the range discussed above. 



<3 < 3 x 10 20 (radiation-dominated) 



or 



a < 10 (matter-dominated) . 
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Anothcr possible source of constraint is the gravitational attraction that 
such a neutrino cloud would create due to its energy density (note that the 
static scalar field contributes here). In section 3 the gravitational acceleration 
due to various clouds was displayed. Nieto, et al (2Cj] have recently discussed an 
anomalous acceleration observed on the Pioneer spacecraft, essentially constant 
from 10 to 50 AU with a valuc of 10~ 9 m/s 2 . Whilc Figure 5 raises the possibility 
of a nearly constant acceleration over a wide range of distances, the magnitude 
would require an average energy, at a density of 2 x 10 15 /cm 3 , of w 50eV, far in 
excess of the values considered here. For this discussion (one generation only), 
this implies that no useful constraints are likely from gravity. On the other 
hand, such considerations add strength to the argument that the range of the 
interaction ought to be of the order of several AU, since the extent of the surface 
is given by that range. Even though the clouds could be much larger than the 
scalar range, such clouds would have a relatively sharp surface and would not 
produce a radial dependence that was gentle enough to appear constant. 

Therefore, assuming the range of the interaction to be of the order of 1 AU 
give or take a few orders of magnitude, since 

1AU w 1.5 x 10 13 cm 

we have 

m s (XAU) « 1.3 x 10~ 18 e^ 
which suggests that we take 

10" 21 < fj, < io- 17 . 

Thus, for a given Kq , we have a range of allowed a, as may be deduced from 
Figure 7. On that figure, the solid vertical line represents the fi appropriate to 
a range of 1AU and m° = 10 eV, the dotted lines on either side representing 
a change of 2 orders of magnitude. The diagonal lines are lines of constant Kq, 
the right hand one for Kq = 4 and the left hand one for Kq = 40, 000. The 
other constraints already discussed appear as horizontal lines representing the 
various upper limits on a. For ju and K in the ranges above, we would obtain 

10~ 42 < a < 10 -31 

Such small values of a trivially satisfy all other operative constraints. 

6 Tritium beta decay 

As remarked in section 5, a scenario which led to a (smaller) cloud of Dirac 
ncutrinos would providc for a more natural explanation of additional counts 
at the end point of a beta decay spectrum than would a scenario involving a 
larger number of Majorana neutrinos. The reason, which is well known pl[ , is 
summarized here. 
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The cndpoint of the spectrum, calculated from the energy release in the 
decay, is not affected by the presence of the background neutrinos. Electrons 
emitted following the absorption of a neutrino will have an energy beyond thc 
endpoint equal to the total energy of the neutrino. Thus, for a cold gas, the 
spectrum will refiect the Fermi-Dirac distribution extending from Eq + m* to 
E a + ep where E is the true endpoint. (One may safely neglect the changc 
in the field energy due to the disappearance of one neutrino.) If the neutrinos 
are Majorana, there will be an equal distribution of anti-neutrinos which will 
block the emission of the lowest energy anti-neutrinos from the normal beta 
decay. Since the distributions are equal, the blocking will exactly balance the 
additional emission, leading to no net increase in the rate. 

In either case (Dirac or Majorana), the change in the spectrum is more 
complex than the simple addition of a spike. The density required to affect the 
experiments produces a Fermi momentum of at least several electron volts, and 
the distribution of additional events must reflect that. It may well be that fitting 
with the correct spectral shape will destroy the feature, reported in [jjj, that a 
fit, equally acceptable to that with a negative mass squared, can be achieved. It 
is also possible that the correct shape will not change or improve the fit. That 
question can only be resolved by fitting actual data with all experimental effects 
included. Note that the analysis described here, as envisag ed by is for one 
neutrino, assumed to be the electron neutrino. We return to this in section 8. 



7 Cloud Dynamics 

To analyze the dynamics of the system of cloud plus star in a galactic envi- 
ronnicnt, and, in particular, to determine if a star stays within the cloud that 
secded its formation, would require a modified Fokker-Planck treatment p2fl . 
While a complete treatment would require the complications of at least three 
generations of neutrinos, we consider only one generation, to illustrate some of 
the issues. 

The primary mechanism for altering a star's trajectory is the gravitational 
scattering between two stars that pass relatively near each other and, to lowest 
order, the clouds simply follow along. Since the cloud-star system is polariz- 
able, there will be an induced dipole-dipole interaction, analogous to atomic 
scattering, which will produce a Van der Waals like interaction. While this may 
alter the specifics of the velocity distribution slightly, it should not have a major 
impact on the issue of the cloud remaining with its star. 

A more serious question involves the interaction between two clouds when 
they touch. If we assume that the cloud contains an energy equivalent to about 
1Mq, uniformly distributed to 10 17 cm., then the gravitational binding energy 
of the star to the cloud is ss 10 39 eV or about 10~ 27 eV/zA According to 
Figure 6 for Kq — 200, the surface contribution to the energy per neutrino is 
« 10" 3 m /N 1/3 . Thus, for the cases of interest here, the surface tension of the 
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clouds overwhelms the gravitational interaction with the stars and the stable 
final configuration would have one star denuded and the other dressed with 
twice as many neutrinos. Binney and Tremaine p2fl present the estimate that 
2 stars, with R w R Q , would actually collide once every 10 19 yr. If, however, 
clouds extend to 10 17 cm. the ratio of geometric cross sections is w 2 x 10 12 , 
so the encounter rate would be about 2 in 10 7 yr., which is relatively fast on 
Galactic timescales. 

Note, however, that the cloud stays with one star or the othcr. The evolu- 
tionary result of such collisions would bc that ncutrino clouds would be found 
only with a fraction of the stars and that that fraction would be smaller in more 
densely populated regions. Furthermore, the simple argument presented above 
takes no account of other neutrino generations. The possibility that the Sun 
has remained with an attendant cloud remains viable. 



8 More Generations 

The preceeding discussion has been confined to the artificial case of one gen- 
eration for pedagogical purposes, but, as the last sections illustrate, it is not 
reasonable to expect a good description of nature without allowing for at least 
three generations of light neutrinos. Necessary as this is, the unknown vacuum 
masses and couplings to the scalar field provide a large number of free param- 
eters. The general arguments setting allowed ranges will not change, but the 
actual limits may depend on the number of participating neutrino species and 
the detailed nature of the coupling of the scalar to those species, consequently 
those limits may change by a few orders of magnitude. On the positive side, 
the existence of three mass eigenstates allows for structures on different scales 
and can avoid the possibilities of apparent inconsistencies described above. 

To illustrate some of the complexity that will arise when all three genera- 
tions of neutrinos are considered, it is instructive to examine the case of two 
generations. In principle, there could be several scalars coupling to the various 
generations with arbitrary strengths, but that problem is too unconstrained to 
be instructive. We shall consider two simple examples, one in which the cou- 
pling is proportional to the vacuum mass and one in which the coupling has the 
same strength to both generations. 

We first consider the case where the coupling of a neutrino to the scalar field 
is proportional to the vacuum mass. Label the two neutrinos by h (for heavy) 
and l (for light). Thcn, if g denotes the coupling constant for h, 

m (0) 

9l = ^M><9 (27) 
m h 

This gives equations for the effective masses 

m* = mf - g<j> (28) 
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and 



mj = mf' - gicj) 
If we now define y h and yi as in section 2, we obtain 

Vh = (i ■ 



m 



(0)' 



(29) 



(30) 



and 



or 



Vi = (l - 



ro, v ro^ 



= y/i 



(31) 



(32) 



That is, there is only one function y for both species. We then define the 
quantities x hl xi, e h and e; in analogous fashion. Also define the ratio 



,(0) 



r = 



,(°) 



(33) 



With these definitions, the differential equation for y becomes 
d 2 y dy K 



where 



F = e Fh x Fh -y 2 ln 



e-F h + x Fh 
V 



+ r* 



eF,x Fl - y 2 \n 



e F , + x Fl 



V 



(34) 



(35) 



The factor of r 4 comes from one power for scaling the energy and three powers 
for scaling thc Fcrmi momentum x Fl . 

The solution of the diffcrcntial cquation is carried out in much the same 
manner as discussed in section 3, with e Fh and e Fl playing the roles of Lagrangc 
multiplicrs fixing the total numbers of heavy and light neutrinos respectively. 

For this system, the density distributions are similar with the species with 
the larger e F extending out to larger z. Of course, the factor r 3 in the density 
and r 4 in both the energy and the driving term of the diffcrcntial equation 
mean that the heavy neutrino dominates the system. For the special case where 
e Fl = e Fh , the densities track exactly and the system is cquivalcnt to one 
generation with K' = (1 + r 4 ) K . The possibility of different values for thc 
e F f allows for the composite distribution to fall more slowly at the surface than 
it would for one generation and, if r is not very different from 1, the gravitational 
force will also fall more slowly. 

An example of this is given in Figurc 8, in which we have chosen the ratio of 
light to heavy vacuum masses, r, to be 1/2 and the remaining parameters to give 
K = 20 for thc light neutrino. We have kept the number of heavy neutrinos, 
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Nh, approximately constant and varied the numbcr of light neutrinos, Ni, from 
to a value slightly in excess of Nh- Since both neutrino densities are involved 
in the driving term given in Eq (35), e Fh must be varied as well as e Fl , hcncc 
the approximate constancy of Nh- For each case we have plotted the individual 
number densities as well as the gravitational forces, calculated as in Figure 5, 
as a function of the scaled radius. 

There is a particularly interesting set of configurations in which e Fh is less 
than the value of y corresponding to the infinite system of light neutrinos with 
a given e Fl . In these cases, which can occur as the Universe expands, outside 
the cloud of heavy neutrinos y approaches this asymptotic value from below. 
These configurations correspond to the formation of clouds of heavy neutrinos 
in the infinite sea of light neutrinos at an epoch in which the light ncutrinos are 
still very relativistic. 

The case in which the scalar couplcs to all generations with the same strength 
allows for more radical changes in the density distributions, and we now examine 
that situation for two generations. Here, 

m* = m (0) - g<j> (36) 

for both generations, i.e. the magnitude of the shift is the same for both. Hence, 
if the shift is significant for the heavy neutrino it may cause the effective mass 
of the light neutrino to change sign. This is like the case of a massless fermion 
gaining mass from the presence of a scalar condensate in field theory; the density 
of heavy neutrinos serves as the source of an external scalar ficld for thc light 
neutrinos (and vice versa, although with little effect). Care must be takcn in 
the assignment of quantum numbers, but the energy, given by 

E = (m* 2 + k 2 ) 1 ' 2 

is well behaved. If y denotes m* h /m h °\ then 

yi = 1 - (1 - y)/r (37) 

ei = {y 2 +x 2 f' 2 (38) 

The solution of the differential equation proceeds as before, with one sub- 
sidiary condition. Since e; > 1 implies that it is energetically favorablc for 
the light neutrino to be at infinity, light ncutrinos will be excluded from such a 
rcgion and x Fl is sct to unless the condition 

- e Fl <yi < e Fl (39) 

is met. For r sufficiently small, one may obtain a solution in which the heavy 
neutrinos occupy a sphere with radius a; from a to some larger radius, b, there 
are no neutrinos and <fi is a linear combination of modificd spherical Bessel 
functions of the third kind of order 0; then light ncutrinos occupy a spherical 
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shell from b to some large radius c. An example of such a distribution is plottcd 
in Figure 9, along with the gravitational force, as a function of the scaled radius. 
Here the vacuum mass ratio was r = 0.1, Kq = 400 for the heavy neutrino 
and JV, « | N h . 

In principle, if the shell of light neutrinos occurs at a large enough radius, 
such a distribution could radically affect the radial dependence of the gravita- 
tional force produced by the neutrinos. In practice, that will require a very 
delicate matching in the solution of the differential equation (troublesome for 
theorists, not for nature) which requires the development of techniqucs like those 
reported is section 2 for the treatment of very large clouds for one generation 
only. 

Sincc this scalar field may well be a contributor to, but not the only con- 
tributor to, the neutrino mass matrix, the real condition is probably more com- 
plicated than eithcr of these examples. The general feature that the heaviest 
mass eigenstate should be most dense and of smallest extent, common in both 
special scenarios, ought, however, to persist. If the dominant component of a 
neutrino cloud in the region of one to several AU from the Sun were one of the 
heavier mass states, say the medium mass eigenstate, with a small component 
to be the interaction eigenstate associated with the electron, it could be possible 
to provide enough energy density to affect satellite behavior without violating 
the bounds from Supernova 1987a or conflicting with Tritium beta decay ex- 
periments. In the latter case, while the effect could extend for several tens of 
electron volts to either side of the end point, a small mbdng angle would provide 
only a small Pauli blocking below and a concomitant small count rate above. 

Whatever the detailed couplings might be, the effective mass of each of the 
mass eigenstates will be affected by the distributions of all three. This could, 
therefore, alter the usual analyses of the MSW effect in the sun and, since the 
neutrino densities associated with supernovae is so very large, could lead to very 
different effects there. Recent analyses of the abundances produced in a neutron 
rich exterior of supernovae [ p3| seems to favor an unusual spectrum of masses, 
which might be achieved, for effective masses, by the scalar field discussed here. 

9 Conclusions 

We have applied the techniques of Quantum Hadrodynamics to the study of a 
system of neutrinos interacting through a light, weakly coupled scalar boson. 
We have shown that, for a wide range of parameters, neutrinos will tend to 
condense into clouds, with dimcnsions the scale of the inverse boson mass. In 
fact, for parameters which cause no conflict with laboratory measurements, such 
clouds could easily be the right size and density to affect experiments on and 
around the earth. 

We have shown that it is likely that any such condensation would have 
occurred before recombination and that the formation of neutrino clouds could 
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form a natural seeding mechanism for the formation of hadronic objects on 
thc scalc of stars. Ncutrino cloud formation, being a phasc change, occurs 
very quickly, so these seeds are available at the earliest possible epoch for star 
formation. 

The extension of this work to more than one neutrino flavor depends on 
the mass hierarchies of the neutrinos and the scalars as well as the details 
of the coupling of each scalar to different generations. For thc casc of two 
generations of ncutrinos and onc scalar field, we have looked at two simple 
choices for the couplings. While different in detail, both generate concentric 
sphcrical (assumed) distributions with the lighter neutrinos (as determined by 
vacuum mass) extending farther out. We would expect this general feature to 
survive for thrcc gcncrations, raising the possibility of the heaviest species being 
essentially within the star, the other two occupying different regions of space 
out to a distance, depending on the detailed history of the system, of a fraction 
of a parsec. 

If thc density of the clcctron component of the neutrinos and antincutrinos 
around the Sun is high enough, there could be observable effects on very sensitive 
experiments such as the study of Tritium beta decay to search for antineutrino 
mass cffects or double beta decay measurements seeking evidence that neutrinos 
are Majorana particles. 

One consequencc of the existence of such an interaction would be that all 
such measurements would have to be interpreted in terms of effective masses, 
rather than the vacuum masses that are relevant to model building. 

Whether terrestrial effects are observed or not, evidence for or against the 
existence of such a scalar interaction is most likely to come from astronomy 
and astrophysics. The implications of the existence of neutrino clouds, with 
respect to both the time scale and the distribution of sizes, should be amenable 
to testing through modelling and observation. The gravitational effects within 
our own Solar system, while subtle, could be observable in vcry high accuracy 
satellite tracking data. Depcnding on thc prccisc modcl for several generations, 
one may be able to observe the modifications of oscillation and propagation in 
the extremely dense neutrino fluxes associated with supernovae. 

Whatever the expcrimcntal outcome of such tests may be, this problem re- 
mains as a fascinating extension of the theoretical techniques of QHD, dcvclopcd 
for the study of atomic nuclei, to vastly different regions of parameter space. 

This work was supported in part by the US Department of Energy, thc 
National Science Foundation and the Australian ARC. Somc of the work was 
carried out during visits to the Institute for Nuclear Theory at the University 
of Washington; their hospitality is gratefully acknowledged. One of us (GJS) 
particularly wants to thank C. Horowitz for several very useful conversations. 
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Figure Captions 



Figure 1. Scaled energies vs. scaled Fermi momentum for infinite matter at 
zero temperature, K = 20. The energies and the Fermi momentum are given 
in units of the neutrino vacuum mass, mi°\ K is defined following Equation 
(7). The dashed curve is the Fermi energy and the solid curve is the average 
total energy per neutrino, which includes the contribution of the static scalar 
ficld, e = e„ + e s . 

Figure 2. Scaled total energy per neutrino vs. scaled Fermi energy for infinite 
matter at zero temperature, K = 20. 

Figure 3. Scaled total energy per neutrino vs. scaled Fermi energy for infi- 
nite matter at zero temperature for several small values of K . For values of 
K > 3.27, the system is self-bound. For 2.67 < K < 3.27 there is a local 
minimum above e = 1, allowing for the possibility of metastable states. 

Figure 4. Number dcnsity p vs. scalcd radius for thrcc diffcrent sized clouds 
for K = 20. In each case ep was chosen, a solution was obtaincd as dcscribed 
in the text and the total numbcr N was obtained by integrating the density. The 
radius is given in units of l/m s ; the number in units of /i~ 3 . For the solid curve, 
eF = -74 and N = 2.881; for the dashed curvc, ep = .76 and N = .851; 
and for the dot-dashcd curvc, ep = -78 and N = .373. 

Figure 5. Gravitational force, in arbitrary units, for the three cases shown in 
Fig. 4. f g at z is given by 4j j Q z e z 1 dz. 

Figure 6. Dependence of the average energy per neutrino and of the radius on 
the number of neutrinos in the cloud, shown for K = 200, on a log-log plot. 
The open squares are the values obtained for the logarithm of the difference 
between e for the cloud and its value for infinite matter; the solid line has slope 
— 1/3. The stars are values of the logarithm of the radius at which the neutrino 
density goes to zero; the dot-dash line has slope 1/3. The filled circles are the 
values of the e^ associated with each calculation, the dotted curve is just a 
guidc for the eye. 

Figure 7. Allowed regions of the /j, - a plane. The horizontal lines indicate 
upper limits on a obtained from the labeled processes as discussed in the text. 
The solid vertical line corresponds to the value of /j for m s ~ 1.3 x 10~ 18 eV, 
appropriate to a range of 1AU, and m' = WeV. The dashed vertical lines cor- 
respond to changing /j by two orders a magnitude each way. The diagonal lines 
represent constant K , the lower for K = 4 and the upper for K = 40, 000. 
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Figure 8. Densities and gravitational force for two neutrino mass eigenstates 
with g oc m. For this example, — = \ and K = 20 for the light neutrino. 
The four cases keep N h approximately constant and vary N. The coupling of 
both densities to the scalar field requires e Fh to vary, hcncc thc approximate 
equality. Case A: N h = .040, N = 0, and e Fh = .75. Case B: N h = .042, 
N = .011, e Fh = .725 and e Fl = .825. Case C: N h = .040, N = .027, 
e Fh = .7 and e Fl = .93. Case D: N h = .041, N = -051,e.F h = -675 and 
e F[ = .975. For each case, the long-dashed line gives the density for the heavy 
neutrino, the dot-dashed line gives the density of the light neutrino and the 
solid line gives the gravitational force due to the total distribution calculated as 
for Fig. 5. While the units of the gravitational force are arbitrary, the relative 
strengths scale properly. 

Figure 9. Densities and gravitational force for two neutrino mass eigenstates 
with constant coupling to thc scalar field. Kq = 400 for the heavy neutrino 

and — = .1 The curves have the same meaning as in Fig. 8. e Fh = .7, 

e Fl = 1 .85, N h = .000339 and N = .00162. 
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